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Statistical thermodynamical and kinetically-limited models are applied to study the origin and 
evolution of space charges and band-bending effects at low angle [001] tilt grain boundaries in 
YBa2Cu3C>7 and the effects of Ca doping upon them. Atomistic simulations, using shell models 
of interatomic forces, are used to calculate the energetics of various relevant point defects. The 
intrinsic space charge profiles at ideal surfaces are calculated for two limits of oxygen contents, i.e. 
YBa2Cu3C>6 and YBa2Cu3C>7. At one limit, Os, the system is an insulator, while at O7, a metal. 
This is analogous to the intrinsic and doping cases of semiconductors. The site selections for doping 
calcium and creating holes are also investigated by calculating the heat of solution. In a continuum 
treatment, the volume of formation of doping calcium at Y-sites is computed. It is then applied to 
study the segregation of calcium ions to grain boundaries in the Y-123 compound. The influences 
of the segregation of calcium ions on space charge profiles are finally studied to provide one guide 
for understanding the improvement of transport properties by doping calcium at grain boundaries 
in Y-123 compound. 



I. INTRODUCTION 



The study of grain boundaries of YBCO is very impor- 
tant, because they act as strong barriers to current flow in 
YBCO This behavior of high temperature supercon- 
ductors is in sharp contrast to metallic low temperature 
superconductors, in which grain boundaries are not only 
transparent to current, but also significantly contribute 
to flux pinning. The formation of such weak links in high 
temperature superconductors is due to the low carrier 
density and strong dependence of T c on hole concentra- 
tion. There are many experiments showing that the crit- 
ical current density through grain boundaries is verysen- 
sitive to the orientation of the adjacent crystallines 0,0- 
Generally, a large misorientation angle will reduce J c sig- 
nificantly. As we know, there are quite rich phenomena 
occurring in grain boundaries: variations in the chemi- 
cal stoichiometry, dissociation of grain boundary disloca- 
tions, impurity segregation, etc. For instance, recently 
Kung et al 3 report that small angle [001] tilt grain 
boundaries in YBCO with (110) planes exhibit partial 
grain boundary dislocations separated by stacking faults. 
The dissociated grain boundary structures have twice the 
number of grain boundary dislocations and shorter inter- 
dislocation core channel widths than Frank's geometry 
rule predicts (in Frank's rule the separation between dis- 
locations of small angle 9 boundaries is proportional to 
fe). This requires further development of the above model 
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of supercurrent flow through arrays of small angle grain 
boundaries. Furthermore, the situation gets more com- 
plicated because grain boundaries can also act as sinks 
and sources for vacancies. In ionic solids this can give 
rise to charging effects. The existence of space charge 
regions near free surfaces in an ionic solid was first pos- 
tulated by Frenkel In an ionic solid, the vacancy 
concentrations at surfaces and grain boundaries are de- 
termined by the individual free energy of formation. The 
cation and anion vacancies usually have different forma- 
tion energies. However, the vacancy concentrations in the 
bulk are determined by the condition of charge neutral- 
ity and usually they are different from those at the grain 
boundaries where strict charge neutrality conditions are 
relaxed. Thus, the vacancy concentration changes as a 
function of distance from the grain boundaries. Usually, 
vacancies in ionic solids are charged defects. Spatial dis- 
tributions of these charged defects lead to a potential 
difference between the bulk and free surfaces and grain 
boundaries 0, IE 13- There are two kinds of space 
charges. Here we define a space charge arising from dif- 
ferences in point defect formation energies as "intrinsic 
space charge" to contrast with those due to segregated 
dopants and grain boundary nonstoichiometries. 

The topic of space charges at interfaces arising from 
point defects has not been investigated so far for 
cuprates. Since holes carry a charge, it is expected that 
a space charge will produce a significant effect on hole 
distributions around grain boundaries for the metallic 
phase of YBa2Cu307_^. Actually the holes can seg- 
regate around charged boundaries even if not metallic. 
These can be localized hole distributions. The generic 
phase diagram of the cuprates shows a wide variety of 
different behaviour at different temperatures and lev- 
els of doping [H EH . At zero doping the cuprates are 
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all insulators, and below a few hundred kelvin they are 
also antiferromagnets (i.e. the electron spins on neigh- 
bouring copper ions point in opposite directions). How- 
ever, when the doping is increased above a critical value 
(about 5%, although this varies from compound to com- 
pound), the antifcrromagnetic state disappears and we 
enter the so-called underdoped region. As the doping 
is increased further, superconducting region is reached. 
For YBa 2 Cu307_(5, the doping is closely related to the 
oxygen content. Here we examine two limits: 8 = 0,1. 
When 8 = 1, this is an insulator, which can be approxi- 
mately treated as an ionic solid. However, when 8 = 0, 
the phase is superconducting, and has mobile holes. The 
existence of holes can add more screening than in the 
insulator case. This is quite similar to the two limits 
of semiconductors: point defect limit and electron/hole 
limit. As the oxygen content changes from O7 to 06, the 
c-axis expands fl2|. The lattice strain has a large effect 
on point defects |13|. Besides, the valence of chain cop- 
per ions decreases from +2 to +1 as the oxygen content 
is reduced. Both factors lead to a dramatic change of 
point defect formation energies in 123 compounds. Con- 
sequently, this will affect the intrinsic space charge. 

It has been known since the early days of YBCO that 
doping with calcium is one way of adding more holes to 
the material. Divalent calcium is roughly the same size as 
trivalent yttrium and so readily replaces it. If a deficiency 
of holes causes grain boundaries to have weak coupling, 
then the critical current might be improved by doping 
calcium into grain boundaries. Note that if calcium is 
doped inside the bulk as well, the bulk T c will be lowered. 
Therefore, doping calcium at the grain boundaries can 
both repair transport properties of grain boundaries and 
keep the high T c value of the bulk. Indeed, Hammerl et al 
[Tij recently developed one beautiful process to increase 
the critical current density above 10 5 (A/ cm 2 ) at temper- 
atures near the boiling point of liquid nitrogen by dop- 
ing Ca into only high angle grain boundaries of YBCO. 
So, the current key issue is to improve transport proper- 
ties near low-angle [001] tilt grain boundaries, which are 
the main t ype s of defects in polycrystallme YBCO wires 
and tapes |l5l [3, 0| . For small angle grain boundaries, 
Gurevich et al 18] argue that the critical current density 
dependence is mostly determined by the decrease of the 
current-carrying cross section by insulating dislocation 
cores and by progressive local suppression of the super- 
conducting order parameter near grain boundaries as 9 
increases. However, space charge effects due to Schottky 
disorder are not considered in their model, although there 
exist clear space charge profiles in cuprate superconduc- 
tors [l^. It is thus desirable to make some theoretical 
inputs for this topic, which can help to provide useful 
guides to find practical ways to improve the transport 
properties of grain boundaries in YBCO. 

In this paper, first we construct a statistical thermody- 
namic model to study the origin and evolution of space 
charges due to the spectrum of point defect formation en- 
ergies as a function of the variation of oxygen content in 



YBa2Cu307_,5. Then we develop one kinetically-limited 
model, including effects of doping and solute segregation, 
to study the process of Ca doping upon low angle tilt 
grain boundaries in YBCO and consequential influences 
on the space charge profile. 



II. SPACE CHARGE PROFILE FOR THE 
INTRINSIC AND DOPING CASES FOR 123 
COMPOUNDS 

The simplest imperfection in a crystal lattice is a lat- 
tice vacancy, which is a missing atom or ion. Electri- 
cally neutral and structure-preserving groups of vacan- 
cies, known as Schottky defects, are required in com- 
pounds. No matter by what means a Schottky defect 
is made, it is necessary to expend a certain amount of 
work per atom taken to the surface. Therefore, the en- 
ergy of the crystal is increased. At a finite tempera- 
ture, the average energy of formation of the group of 
Schottky defects determines the concentration of such 
point defects at equilibrium. For a multicomponent crys- 
tal, the difference of the energy of formation of various 
defects determines the nature of space charge near in- 
terfaces 5]. We calculated the point defect formation 
energies for both YBa^u^Ofiand YBa 2 Cu 3 07 by the 
Mott-Littleton approach 120, [21j which has been imple- 
mented in GULP codes |22j. The pair potentials of the 
shell model used in this study were developed by Baet- 
zold for Y-123 for studying ionic and electronic (polaron) 
defects |23|. The calculated point defect formation ener- 
gies are collected in Table ©. The data for YBa 2 Cu 3 7 
agree nicely with the previous results [23| using HADES 
III codes M- It can be seen from Table (I) that there 
is a wide spectrum in the calculated formation energies 
of individual vacancies in the Schottky group for both 
YBa 2 Cu306 and YBa 2 Cu307. Thus, in principle there 
should be space charges near free surfaces at equilibrium 
or quenched in from high temperatures where ionic mo- 
bility is high enough to maintain equilibrium. A similar 
phenomenon is expected to occur near internal interfaces 
such as grain boundaries, but the situation is more com- 
plex than for free surfaces for at least two reasons. First 
of all, the spectrum of dopant formation energies may be 
different for the process of removing an ion from a site 
in the bulk crystal and placing it at a site in the grain 
boundary. Of course, at equilibrium the chemical poten- 
tial of defects is the same throughout the crystal, but the 
local variations in concentration will be affected by dif- 
ferent energies for excess species (i.e. the ions removed 
to form the vacancies and placed at surface or boundary 
sites) in different types of interfaces. The second compli- 
cating factor is that the local stoichiometry in the struc- 
tural units (or individual dislocation cores) making up 
the boundary may be significantly different from that of 
the bulk in complex ionic crystals. Advanced transmis- 
sion electron microscopy has been used recently to show 
that this is the case for grain boundaries in SrTi0 3 [2{| 
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TABLE I: Point defect formation energy (eV) in YBa2Cu3C>6 
and YBa2Cu307. This defect formation energy refers to re- 
moving one ion from its lattice site and putting it at ledge- 
corner site on the surface of the crystal. 





YBa 2 Cu30 6 


YBa 2 Cu 3 07 


E(Y) 


14.42 


12.71 


E(Ba) 


4.06 


4.25 


E(Cu-chain) 


0.96 


4.23 


E(Cu-plane) 


6.92 


5.68 


E(O-plane) 


3.47 


3.56 


E(O-chain) 




4.22 


E(O-apex) 


6.85 


4.48 



and for twin boundaries in BaTiC>3 |26j . It is not entirely 
clear how such local chemistry at boundaries in complex 
oxides will affect the intrinsic space charges arising from 
the spectrum of defect formation energies in equilibrium 
groupings, e.g. Schottky defects. Therefore, we will re- 



strict our detailed calculations of intrinsic space charge to 
free surfaces, but we expect the behavior at grain bound- 
aries to be at least qualitatively similar. 



A. Intrinsic Space Charge at Surfaces 

1. Insulator Case: YBaiCuzOe, 

For the sake of simplicity, we assume the crystal has 
free surfaces at x=0 and x=2L and is of infinite extent in 
the y and z directions in studying intrinsic space charge 
problems, i.e. only containing free surfaces. As discussed 
above, the case of grain boundaries is more complex, but 
the calculations for free surfaces will provide some qual- 
itative guide to understanding of intrinsic space charges, 
i.e. those due to equilibrium Schottky defects at other 
interfaces too. The free energy per unit area of a slab-like 
disordered crystal of the thickness 2L is given by 
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r 



where is the electrostatic potential, p(x) is the lo- 

cal charge density, and n Y , n Ba , n C u c , n Cup , n 0p , and n 0a 



yttrium, barium, chain copper, plane copper, plane oxy- 
gen, and apex oxygen sites, respectively. Note that the 



stand for the densities (per unit volume) of vacancies at configurational entropy density S con f is calculated as 
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k n Y n Ba n Cuc 



2N 4N 2N 

<5n Cu In h <m 0p In h <m 0a In , (2) 

n Cup n 0p n Qa 



where N is the number of formula units per unit volume. 
This arises because we have neglected effects of clustering 
and association (e.g. the formation of di-vacancies etc.) 
among the defects. The charge density p(x) is given by 



The electrostatic potential $(x) is obtained from Pois- 
son's equation, which can be written Using equations 10 
HJ) as 



p(x) =e(2no p +2no a -3n Y -2nBa-ncu c -2ncu p )- (3) 

From SF — 0, the spatially- varying equilibrium con- 
centrations are 
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TABLE II: Bulk Potential and Debye Length of YBa 2 Cu 3 6 . 



T (K) 


kT 


k (m) 


350 


-27.42 


44.34 


400 


-23.81 


0.93 


600 


-15.61 


1.83xl0~ 4 


1000 


-9.28 


2.49xl0" 7 



where e is the static relative permittivity (e = 5.21), 
which is obtained by atomistic simulations using our 
shell model I^tI with pair potential parameters from 
Baetzold|23|. We define the electrostatic potential to be 
zero at the surface, and the bulk potential $oo is de- 
fined as the potential for which p = 0. Therefore, letting 
the right-hand-side (RHS) of Eq. (J5J equal zero, the bulk 
potential can be obtained, and the defects' concen- 
trations in bulk can be calculated. We also define the 



Debye length k [E IE as 



where nj is the density of vacancies on sites of type i, and 
qi is the charge of the vacancy. 

If we scale the original length r by the Debye length k, 
and define one new variable z as follows, 

x 

s = — , 

K 

kT ' { ] 

equation J^J can be rewritten as 
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Boundary conditions for the slab of thickness 2L are given 
by 



s = 0, z = - 
L e($ L 



e($oo) 
kT ' 



kT 



(9) 
(10) 



where is the potential at x=L. In this study, we as- 
sume the thickness approaches infinity, which means the 
grain boundary is qualitatively modeled by two surfaces 
"back to back". Under this condition approaches 
<&oo. We solved the Eq.© by conventional numerical 
integration techniques for several temperatures. The re- 
sults are plotted in Fig^ Since there is no oxygen at 
chain sites for YBa2Cu30g, the chain copper is compar- 



atively loosely bound inside crystal. Based on the point 
defect energy calculations, it turns out that it costs the 
least energy to make such a copper vacancy (see Table 
(P) ) . Because the free energies of formation of the anion 
and cation vacancies differ, there exists a charged sur- 
face and a region of space charge of the opposite sign 
adjacent to the surface. For YBa2Cu3C>6, this leads to a 
positively charged surface and a region of negative space 
charge beneath the surface. Presumably this also would 
be the case for dislocations and grain boundaries. 

Based on the point defect energies in YBa 2 Cu 3 06 (Ta- 
bic (Qi), the Debye length k and bulk potential $oo are 
tabulated in Table {nj for various temperatures. As the 
temperature rises, the concentrations of point defects in- 
crease exponentially such that the Debye length becomes 
reduced at higher temperature. Next, let's consider the 
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FIG. 1: Free surface space charge in YBa2Cu30e. $00 refers 
to the bulk potential caused by space charge. S is the reduced 
length, scaled by the Debye length. 



trend of bulk potential versus temperature. In the RHS 
of Equation JSJ , it is seen that at lower temperature there 
is a bigger difference among exp(— ) . To balance this, 
a larger absolute value of the bulk potential relative to 
the surface is needed. Since this potential is positive, 
this expels holes from the surface. In other words, the 
positively charged surface can cause a serious depletion 
of the hole content there. 



2. Metallic Case: YBa 2 Cu 3 07 

Compared with YBa2Cu306, the defect energy "spec- 
trum" as calculated with the shell model with pair po- 
tential parameters from [23j is different for YBa2Cu3C>7 
(see Table ([}). The similar comparison has also been 
discussed by Baetzold in studying point defects in 
YBa 2 Cu307 and YBa 2 Cu306.5 |23H28||. Since the plane 
copper is coordinated with four oxygens, it is quite sta- 
ble in the sense that more energy is needed to make a 
vacancy at the copper site in Cu02 planes. Interestingly, 



our atomistic simulations show that the formation en- 
ergy of oxygen vacancies in the CuO-2 plane is the lowest 
one in YBa2Cu30y. Another very important difference is 
the existence of screening by mobile holes in YBa 2 Cu307, 
which is absent in YBa2Cu306 because it is an insula- 
tor. In our atomistic simulations of the formation energy 
of point defects we ignore the explicit screening of de- 
fects by mobile holes. However, the screening by holes 
is included in the space charge calculations, assuming 
holes are localized on ion sites (but mobile by hopping 
between ions). Considering the temperature of fabrica- 
tion of YBa2Cu307 is much higher than T c , YBa2Cus07 
is in fact in a normal state. The Ioffe-Refel resistivity 
of YBa2Cu307 at normal state is very large, about 0.1 
milcm, due to the small carrier concentration [2^. l8fl | . 
Emery and Kivelson suggested this type of exotic "bad 
metal" can be still regarded as a quasiparticle insulator 
|3lj . These observations encourage us to apply the clas- 
sical approach to study the effects of band-bending on 
hole's segregation near surfaces at fairly high tempera- 
ture (above 600K), which YBa2Cu3C>7 remains in metal- 
lic normal state. 

Following a procedure similar to that used for 
YBa2Cu306, but now including chain oxygen vacancies 
and localized (but mobile) holes, from SF = the equi- 
librium concentrations are 
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= 2Nexp(- 


ncu c 
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= 2Nexp(- 
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= 4Nexp(- 
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= 2Nexp(- 


n o c 


= Nexp(- 



(11) 



E, 



kT ^ ' 
+2e# \ 



Eo e +L$ 
kT 



where no c stands for the density (per unit volume) of va- 
cancies at chain oxygen sites. Since mobile holes are 
present in YBa2Cu307, the charge density is given by 



p(x) = e(2no p + 2n 0a + 2n 0c - 3n Y - 2n Ba - n CUc - 2n CUp + Dhole) < (12) 

I 



To calculate the bulk potential $00, we assume that 
holes are bound to oxygen ions in the Cu02 plane such 
that the effective valence of oxygen changes from —2 to 
— 1.75. Within this approximation, we compute the in- 
dividual defect energy by the Mott-Littleton approach. 
Then by applying the charge neutrality constraint for the 
bulk, the bulk potential <&oo is obtained. The resulting 
values of bulk potential and Debye length are tabulated 
in Table 



Notice that when including the screening effect caused 
by mobile holes, the bulk potential is invariant. However, 
the Debye length changes dramatically because of the 
hole concentration. At the optimally doped case, there is 
around 0.2 hole per Cu02 [13 ■ Therefore, the total hole 
content is about 0.4 per unit cell of YBa2Cu307. We 
assume that holes are localized, but, can hop from site 
to site. Hence, the hole content at the position of x with 
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TABLE III: Bulk Potential and Debye Length of YBa 2 Cu 3 7 



T(K) 



k (m) 



600 3.55 l^lO-" 11- 
800 2.76 1.82xl0 -11 
1000 2.16 2.04xl0" n 



where N-fo is the density of holes at surfaces, determined 
by the following condition 



l).4xN Y / (1A / dxN-f xexp ^°°) . (14) 



the electrostatic $(x) inside the slab is given as 



n holc (x)=N-f xexp $oo) (13) 



where dA = dy * dz is the area. After the same process 
as described in the previous section, the Poisson equation 
yields 



-) exp(— 2z) 
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(15) 




s ( X / K ) 

FIG. 2: Free surface space charge in YBa2Cu307. 3?<x> refers 
to the bulk potential caused by space charge. S is the reduced 
length, scaled by the Debye length. 



In YBa2Cu3C>7, based on the shell model calculations 
of point defect energies, it turns out that it is much 
easier to create oxygen vacancies in the CuC>2 plane 
than any other type of vacancy (See Table (JU). The 
space charge profile is obtained by numerical solutions 
of Eci. (|15fl at several temperatures, which are plotted in 
FigEl In our model, in contrast to YBa2Cu306 (insula- 
tor), YBa2Cu307 is a p-type superconductor, containing 
mobile holes. Assume that the total hole content is 0.4 



per unit cell, this makes the dominant contribution to 
the total charge density such that the total charge den- 
sity is effectively constant over the temperatures range of 
interest. As for the bulk potential relative to the surface 
calculated by the procedure described in the previous 
section its value is determined by differences among the 
cxp(— terms. So, at higher temperature, a much 
smaller bulk potential is needed to keep the bulk neu- 
tral. As the total charge density is effectively constant, 
the Debye length is proportional to the square root of 
the temperature [see Eq. (6)]. It is intrinsically different 
from the corresponding trend for YBa2Cu3C>6. Further- 
more, note that the Debye length of O7 is far shorter 
than that of 06, the mobile hole provides strong screen- 
ing of the charge at grain boundaries. We also computed 
the Thomas-Fermi screening length Atf by the following 
relation [32l |: 

h , £ n 1 / 7T ,1 , „ 

A TF = -( )2 - 6 (16) 

eMmhole 3 x (0.4 x N) ; ' v ' 

where mhoic is the effective mass of the hole, N is the 
number of YBa2Cii307 unit cell per volume, and we as- 
sume the total hole density equals to 0.4 x N for optimal 
doping. The calculated Atf is around 1 A by taking 
nihoio = rn c i c . It is interesting to note that both lengths, 
the Thomas-Fermi screening length Atf and the Debye 
length k, are shorter that an interionic spacing. Because 
holes are localized, the effective mobile hole density is 
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expected to be only small amount of the total density. 
In fact, the k and Atf are the lower limit of the effective 
screening length. The very short screening length can 
help to protect superconducting properties inside bulk 
from the influence of grain boundaries if the hole density 
of grain boundaries corresponds to YBa2Cu307. How- 
ever, experiments have shown that the oxygen content at 
grain boundaries is often less than that inside the bulk [2] . 
If oxygen content is less than around O6.45, from the cal- 
culation in the previous section it can be concluded that 
the Debye length will increase dramatically. It means 
that a large portion of the bulk near the boundary will 
be affected by the internal electric field caused the space 
charge around grain boundaries. This may be one rea- 
son for the vicinity of some grain boundaries not being a 
good superconductor. 

B. Space charge of calcium-doped YBa2Cu3C>7 

After the discovery of YBa2Cu307, many types of ions 
have been chosen to dope into YBa2Cu30y in search of 
new superconductors and to improve properties. Dop- 
ing with Ca into [001] tilt grain boundaries of YBCO 
has been shown to increase the critical current density 
above 10 5 (A/ cm 2 ) 14]. Our goal here is to study the 
segregation of calcium around small angle tilt bound- 



The first case, Eq. l|17|l . is homovalent substitution. 
Simply, barium is replaced by calcium. However, yt- 
trium is replaced in the last two reactions fEas. (|18|l and 
H19J0 . These two reactions are aliovalent substitutions, 
in which it is very important to identify the mechanism 
of charge compensation. In reaction (|18fl . a hole is cre- 
ated in the CuO chain to make the whole system charge 
neutral. Assuming the hole is trapped at copper site, 
the valence of one copper at chain site increases by one 
approximately. In reaction i|19|) . we assume that the ex- 
tra hole is in the Cu02 plane, such that the valence of 
one copper in the plane increases by one. All the cal- 
culations are done by assuming the infinite dilution of 
dopants. When the valence of copper changes, the corre- 



aries of Y-123 in order to aid in the search for a more 
efficient approach for doping Ca. First, doping calcium 
in bulk is studied by shell model calculations to iden- 
tify the energetically favorable reaction for doping. Sec- 
ondly, the volume of solution of a Ca dopant atom is 
calculated. Thirdly, the segregration of calcium is repre- 
sented in terms of the segregation energy as function of 
angle around the core of edge dislocations in tilt bound- 
aries, which can clarify the role of the elastic strain field. 
Finally, the space charge profile is studied by including 
extra holes due to dopants and segregation of dopants 
due to the strain field. 



1. Site Preference 

When Ca replaces Y, it is accompanied by a hole for 
the sake of charge neutrality. However, when substituting 
Ca for Ba, the charge remains balanced without chang- 
ing the total hole content. Thus, it is important to know 
the circumstances when Ca doping yields Ca on Y or Ca 
on Ba sites. Here, we mainly focus on the dissolution 
of divalent Ca 2+ , which can be represented by the fol- 
lowing defect reactions, depending on which sublattice 
Ca 2+ chooses and what kind of charge compensation is 
involved. The parameter m is chosen to characterize the 
dopant concentration in the YBCO. 



(17) 



(18) 



(19) 



sponding interatomic pair potentials are scaled based on 
a simple radius scheme |l3j . GULP was used to compute 
the defect energies. The calculated energies of solution 
(per calcium) are listed in Table Previous studies 

on bipolaron binding energies have shown that the hole is 
more likely to localized on copper sites than oxygen sites 
in the Cu02 planes |2fl| . In our studies, the hole content 
is tuned by the dopant, calcium. The calculations sug- 
gest that it is most energetically favorable to substitute a 
calcium ion at an yttrium site, together with the creation 
one hole in the Cu02 plane. This agrees very well with 
experimental measurements jl4j . 



I 

CaO + mYBa 2 Cu 3 7 

-> mY(Ba^-i Ca 1 hCu 3 Q 7 + BaO\ 

2CaO + 2mYBa 2 Cu 3 7 + \p 2 

2m(Y m -iCa ^)(BaO)?(Cu 3+ Q) ^ (CuO) a? -i (CuO^ + Y 2 Q 3 ; 

2CaO + 2mYBa 2 Cu 3 7 + ^02 

-> 2m{Yr^Caj_){BaO) 2 {CuO){Cv 3+ 2 )^_{Cu0 2 ) i ^i + Y 2 3 ; 

I 
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TABLE IV: Calculated Energies of Solution of Calcium at 
Infinite Dilution in YBa2Cu307 (in eV per Ca 2+ ) 

Reaction l)l7fl 8.14 
Reaction JTSJ 2.99 
Reaction JTHJ 1.02 



2. Volume of Solution for doping AE(2+) at Y(3+) Site 

Doping Ca in grain boundaries of YBCO can improve 
the transport properties significantly since more holes are 
created near grain boundaries. The strong elastic strain 
field around grain boundaries has an important role in 
the segregation of dopants. The disruption of the lat- 
tice by a point defect, such as a dopant atom like Ca + 
substituting for Y 3+ in YBCO, gives rise to forces on 
the surrounding lattice atoms which, in turn, give rise 
to a displacement field in the lattice around the defect. 
The interaction between the defect displacement field and 
the stress fields of other defects, such as dislocations, 
or applied stress fields, alter the energy of the defect 
[33L 134 and cause defect segregation. The effects 
of the forces on the lattice due to the point defect can 
be described by means of force multipoles |33| , which to- 
gether with the Green's tensor function of the matrix of 
the perfect crystal, can be used to calculate the resulting 
displacement field caused by the defect in the crystal. In 
this section, we use numerical shell model calculations 
with GULP to estimate the strength of the elastic force 
dipole for the point defects (dopants). Then in the next 
section we apply anisotropic continuum elasticity to de- 
scribe the stress field due to an edge dislocation; finally 
we use continuum elasticity theory to compute the elastic 
interaction between the dopant atom and the dislocation. 
Such methods have been shown to give good agreement 
with direct numerical calculations of the interaction be- 
tween external pressure and vacancies in simple crystals 
such as van der Waals bonded rare gas crystals (see ap- 
pendix 3 in reference [3^|L 

To study the segregation of calcium around tilt grain 
boundaries, it is necessary to have the volume of solu- 
tion to compute the contribution to the binding energy 
of Ca 2+ at grain boundaries due to interactions with the 
strain field of edge dislocations in the boundaries. In the 
force-dipole approximation, the difference between the 
energy to form the defect in a crystal containing a strain 
field ejj (either from an applied stress or from the stress 
field of another defect) and the energy in the strain-free 
crystal is given as 



Eint — —ttijCij, (20) 

where the Einstein summation convention for repeated 
indices is used and the coefficients a^ characterize the 
strength of the defect forces on the lattice. The effective 



body force, f, which the defect exerts on the lattice, is: 

/i = -°«^(0, ( 21 ) 

where S(r) is the Dirac delta function. These three 
force components represent three orthogonal "double 
forces without moment" and are sometimes called "force 
dipoles" . The values of ay , which have the unit of energy, 
must be obtained from experiment or numerical simula- 
tions. For a defect with spherical symmetry in an elasti- 
cally isotropic matrix: 

«ij = aSij, (22) 
where Sij is the Kronecker 5 (=1, if i=j; =0, if i^j), and 

a = BQ, (23) 

where is the volume change produced by the formation 
of the defect in a finite crystal and B is the bulk modulus 
[35| . Thus from equation (|20|) we obtain 

E mt = -fln(en + e 22 + £33) = PSl, (24) 

where P is the hydrostatic component of pressure at the 
defect, produced by external forces or other defects. 

In the case of a defect with orthorhombic symmetry, 
equation lj2TJ|) yields 

Eint = -[fflnen + a 2 2£22 + 033633]. (25) 

Clearly, it is necessary to evaluate Ej n t for three uniaxial 
strain cases. For strain along a-axis (en ^ 0; £22 = £33 = 
0), we have 



£11 

where Ej n t is evaluated from numerical simulation by 
GULP. Using generalized Hooke's law in the form of elas- 
tic compliances, Sy , we obtain by choosing the crystal 
axes as principle ones 

£ll = SlifTii + S12O-22 + Sl3ff33, (27) 
£22 = Sl20"ll + •S'220'22 + <S23C33, 
£33 = <5'l3 cr ll + <S'23 cr 22 + S33O33. 

Hence, the components of volume of formation of dopants 
are given as 

^11 = Siicin + Si 2 a 2 2 + Si 3 a 33 , (28) 

^22 = 5*12011 + 522022 + 523033, 
^33 = S13O11 + 523022 + 533O.33. 
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the basal plane isotropy, the stress field is given by 




1.1 1.2 
Radius of A E (A) 



FIG. 3: Components of Volume of Formation for doping 
AE 2+ at a Y 3+ site of YBa 2 Cu 3 7 vs. Radius of AE. 



The calculated components of the volume of solu- 
tion for various divalent ions (AE 2+ ) at the Y 3+ site of 
YBa2Cu307 are plotted in FigEI Considering the trend 
across the series of AE 2+ , the ion size plays the important 
role in the stress effects caused by doping AE 2+ at Y 3+ . 
The point defect for doping AE 2+ at Y 3+ site carries one 
negative charge. This strongly repels the surrounding 
oxygen ions to make the sign of the volume of solution 
positive for AE ions whose radius is similar or larger than 
that of Y 3+ (1.02 A). Note that the short range repulsion 
between AE and surrounding oxygen ions becomes larger 
with increasing AE ion radius, it is expected that the vol- 
ume of solution increases with AE ion radius. However, 
for an AE ion with a larger radius, the polarization of the 
lattice around the solute which helps to lower defect's en- 
ergy by polarizing electron "clouds" also depends on the 
distortion of the lattice by the AE ion and thus con- 
tributes to the volume of formation. Therefore, because 
of these competing effects the volume of formation ends 
up with a "satuation-like" curve. 



3. Dislocation Stress Fields in Y-123 

Grain boundaries of importance in coated conductors 
are low-angle [001] tilt boundaries. Thus they can be 
considered as an array of isolated edge dislocations. We 
observe that the elasticity in YBa2Cu30y is not far differ- 
ent between the a- and b-axis ( (C xx -C yy ) /C yy ~ 0.082), 
so this orthorhombic crystal can be approximated by 
a hexagonal one, which is motivated by the existence 
of closed-form expressions for the stress fields in the 
case of anisotropic elasticity in hexagonal crystals. Here 
we average elastic constants of orthorhombic symmetry 
into those of permitted by hexagonal symmetry by tak- 
ing arithmetical average of corresponding elastic moduli. 
The geometry of the dislocations in [001] tilt boundaries 
in Y-123 is like prismatic slip in hexagonal crystals. The 
Burger's vector lies in basal plane < 1210 >), while 
dislocation line runs along the hexad axis. The prismatic 
planes are < 1010 >. Following Ref. ,36], on account of 



'mi 



K x -b-y (3x 2 + y 2 ) 
2n-(x 2 + y 2 ) 2 ' 
= -Ki-b-y (x 2 - y 2 ) 
2tt • (x 2 + y 2 ) 2 ' 



(29) 

(30) 
(31) 



where s is the elastic compliance modulus, and K\ = 
tt5 ^ — j-r. (Note: The formula of Ki on page 192 of 

Ref. [3(| has one misprint.) 

Small angle (usually less than 5°) grain boundaries can 
be considered as an array of isolated dislocations. The 
separation between dislocations is around d = |. If we 
assume that a superposition approximation can be ap- 
plied here, the total stress field can be written as 



total 



E « 

=-2.. ..+2 



i(y + nd). 



(32) 



When doping calcium at yttrium sites, the non-zero 
tensional elements of the "volume of solution" will in- 
teract with the stress field. Neglecting the effect of the 
solute atom on the local elastic modulii, the segregation 
energy of calcium can be written as 



E. 



i— x : y ,z 



.total 



i2« 



(33) 



In the previous section, we calculated the energy for 
doping calcium at yttrium sites inside a stress free bulk 
Y-123. Here we computed the segregation energy for 
calcium due to the stress field around very small angle 
(usually less than 5°) tilt boundaries, which is plotted 
in FigQ] Since our region of interest is at a small dis- 
tance from the dislocation cores, there exist possible er- 
rors in using linear elasticity theory here and our results 
are probably overestimated. However, the relative energy 
differences are more important than absolute values. In 
Fig0] the calculations show that calcium segregation is 
expected to the tensile stress region due to the positive 
sign of volume of formation. This will affect the dis- 
tribution of mobile hole around grain boundaries. We 
will study this further in the following section for general 
grain boundaries. 



4- Effects on space charge due to segregation of calcium 
near grain boundaries 

It has been noted [33] that edge dislocations could act 
as sources and sinks for vacancies such that space charges 
would appear also in the regions near dislocation cores. 
The modeling of space charges around dislocation cores 
is more difficult than our treatment of the intrinsic space 
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FIG. 4: Elastic components of the segregation energy of cal- 
cium around the center of an edge dislocations in small an- 
gle tilt grain boundaries YBa2Cu307. $ is the angle starting 
from (0,b), rotating counterclockwise around the core of small 
angle (usually less than 5°) tilt grain boundaries for a distance 
of one Burgers vectors from the center. 



charges at free surfaces from the difference in formation 
energy of equilibrium point defects. If for no other reason 
than that different in the energy required for creation of 
point defects by removing ions from the bulk and placing 
them in the dislocation core, it is not likely to be quite 
the same as that for free surfaces. However there is also 
likely to be an alteration of the space charge by dopants 
and their accompanying holes if the dopants tend to seg- 
regate to the dislocations. In the previous section (II 
B.3) it was shown that there is a strong attraction of the 
AE 2+ dopants, i.e. Ca, Sr, and Ba, to the tensile regions 
of isolated edge dislocations in low angle grain boundaries 
and a repulsion from the compressive regions. We will il- 
lustrate some of the effects of such dopant segregation 
on space charges by an idealized planar grain boundary 
which consists of joining together two free surfaces with 
their associated point defect space charges, together with 
Ca 2+ dopants and accompanying holes which have seg- 
regated in the stress fields of the boundary. We assume 
that the dopant atoms are attracted to or repelled from 
region along the boundary by the elastic fields which de- 
cay going away from the boundary. The results of this 
calculation will illustrate the interplay between dopant 
segregation and point defect induced by space charges. 

From previous thermodynamic calculations we know 
that hole density changes with the content of doped cal- 
cium. Consequently, the Debye length varies with Ca 
concentration. Note that the defect created by replacing 
Y 3+ with Ca 2+ carries one negative charge. Assuming 
the total number of Ca 2+ is conserved, after solving the 
distribution of charged particles in an electrostatic field 



caused by space charges, the defect density nc a at the 
position of x with the electrostatic ^(x) away from the 
boundary plane is given as 



n Ca = N ■ f Ca * exp 



kT 



(34) 



where N ■ fc a is the density of calciums at surfaces de- 
termined by the following condition 



^ ' ^doping ' * 



y~dA y'dxN-f Ca -exp( 



kT 



-). (35) 



where n^ ping is the total concentration of doped calci- 
ums. Since each Ca 2+ is accompanied by one mobile hole, 
the density of holes near surfaces (N • /o) is determined 
by the following condition (assuming the total number of 
holes are conserved) 



(0.4+n^ ping )NxV- / dA 



dxN • fg • exp 



-e($-joo) 
kT 



Note that we extend our studies to general grain 
boundaries here based on understanding divalent 
dopants' interaction with stress field of small angle tilt 
boundaries. As there can be a significant strain field 
in the grain boundaries region, originating from size mis- 
match between solute and matrix atoms, dislocation, etc, 
a strain energy term should be included in the studies of 
Ca-doping because the strain field due to dislocations at 
grain boundaries can substantially affect the segregation 
of Ca. We ignore the elastic interactions with the bound- 
aries of the charged vacancies, and approximate the effect 
of the stress field of the boundary by an elastic term of 
the calcium by the same form used in Ref. @: 



= U (l-(^) n ),s- K 



Uclastic - v i m ■ v . 

d 

= 0, s • k > 2d 



< 2d; 



(37) 



where Uo is the binding energy near the dislocation core. 
We choose n=2, and d is about two lattice parameters 
along b-axis. This is consistent with the periodicity of 
the structural units proposed in the coincidence models 
of grain boundary structures |38| . Calculations in the 
previous section indicated that the volume of solution 
for doping calcium at yttrium sites is positive. Hence, 
the positive sign of Uo means in the compression region 
which repels Ca, while the negative sign of Uo means in 
the tensile region which attracts Ca. The charge density 
now is given by 



J 



p(x) = c(2n Qp + 2n Qa + 2n 0c - 3n Y - 2n Ba - n CUc - 2n CUp + n ho ic - n Ca ). 



(38) 
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Following the same process as previous section, the Poisson equation can be re-written as 



-) exp(— 2z) 



(ds 2 ) - 2 [2X4XeXp( kT 

+ 2 x 2 x exp( ) exp(— 2z) + 2 x 1 x exp( ) exp(— 2z) 



kT 



kT 



3 x 1 x cxp( ) exp(3z) — 2 x 2 x exp( ) exp(2z) 



kT 



kT 



2 x 1 x exp( ^— — ) exp(2z) - 2 x 2 x exp( £ — — ) exp(2z) 



kT 

p / \ r t U elastic \i 

+ /o x cxp(-z) - / Ca x exp(z 



kT 



(39) 



When calcium is doped in YBCO, the process is of- 
ten designed so that calcium diffuses into grain bound- 
aries first. However, it is difficult to estimate the amount 
of calcium at the grain boundaries during the deposi- 
tion process. In this section, we assume the hole den- 
sity and the amount of calcium near grain boundaries 
are 0.6 and 0.08 per unit volume, respectively. Fur- 
thermore, the elastic binding energy of calcium is cho- 
sen to include both compression and tension regions: 
U = (-6,-4, -2, 0,2,4, 6)kT. This is shown in Fig.5 
for an idealized boundary in YBa2Cu3C>7. As discussed 
in previous section, the grain boundaries in O7 may be 
negatively charged because of the intrinsic space charges 
due to point defects. It should be also stressed that cal- 
cium is doped at yttrium sites, such that this site has 
one negative charge accompanied by making holes at 
CuC>2 planes. Therefore, the screening length becomes 
even shorter than that of pure YBa2Cu307. If the cal- 
cium is pushed away from grain boundaries, the space 
charge remains almost the same (somehow little depen- 
dent on repulsive energy). However, on the other hand, if 
the calcium is attracted by grain boundaries, extra nega- 
tive charge will accumulate around grain boundaries such 
that the space charge's profile is significantly adjusted. 
This can strongly perturb the distribution of mobile hole 
density around grain boundaries. The stress field around 
grain boundaries is also so strong (around 10 2 GPa based 
on our calculated elastic properties of YBCO that it sig- 
nificantly affects the surrounding oxygen content and or- 
dering (l3j). If oxygen content is close to 06, the grain 
boundaries are possibly positively charged, and the large 
Debye length leads to a larger portion of the bulk being 
exposed into the influence of grain boundaries. Under 
such circumstances, doping with calcium not only neu- 
tralizes (by the intrinsic negative charge of Ca^ 1 ) the 
charged boundaries, but also increases the density of mo- 
bile holes to screen the charged boundaries. Therefore, 
the bad low angle tilt boundaries are "cured" by doping 
calcium. 




— ■— u = 





— .— u = 


-2kT 


U = 


-4kT 


— T— U = 


-6kT 




2kT 




4kT 


u = 


6kT 



FIG. 5: Space charge of doping calcium at idealized 
YBa2Cu307 grain boundaries at 723°C. S is the reduced 
length, scaled by the Debye length. U is energy due to strain 
field. Negative sign means attractive interaction. 



III. CONCLUSIONS 

We have systematically studied the space charge pro- 
file for YBa2Cri306, YBa2Cu307, and Ca-doped cases. 
Intrinsic space charge is strongly dependent on the 
oxygen content of YBa2Cu307_a, which is determined 
by the annealing temperature and oxygen partial pres- 
sure. For YBa 2 Cu 3 6 the chain copper is loosely bound 
with neighboring two oxygens, therefore it is compara- 
tively easier to create vacancies at copper chain sites for 
06- Hence, the grain boundaries are possibly positively 
charged. On the other hand, the calculations show that 
it requires less energy to make vacancies at oxygen sites 
for O7. Consequently, it was seen that there exists a dra- 
matic change in the nature of the space charge in going 
from 06 to O7. It is interesting to note that the origin 
of the screening in going from 06 to O7 is completely 
different. For 06 it is the thermal activated vacancies, 
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while for O7 it is the mobile holes. This is very simi- 
lar to the two screening mechanisms in semiconductors: 
one is intrinsic, the other is free carriers due to doped 
impurities. In addition, We have carried out studies of 
the segregation of calcium at tilt grain boundaries in the 
Y-123 compound. By comparing calculated heat of solu- 
tion for various doping mechanisms, we have found that 
the most energetically favorable defect reaction is that of 
doping calcium at yttrium sites accompanied by creating 
one hole in the CuC>2 plane per Ca ion. This conclusion 
is valid for the reaction in the bulk, which agrees with 
previous experimental measurements |l4j . To study the 
segregation of calcium at tilt boundaries, the volume of 
formation is needed in order to estimate the binding en- 
ergy to the elastic strain field which arises from the grain 
boundary dislocations. We computed this for the doping 
series of alkaline earths at yittrium sites. After including 
the strain field contribution, the segregation energy of 
calcium at small angle tilt boundaries is obtained. The 
theoretical results suggest that more calcium is expected 
in the tensile regions around tilt grain boundaries. For 
calcium-doped YBa2Cu30y_5 case, our calculations in- 
dicate that Y 3+ ions are replaced by Ca 2+ ions, accom- 
panied by creating holes in CuC>2 plane. The increase 
of hole content strongly enhances the screening effect by 
mobile holes. The strong screening effect (much shorter 
screening length) helps to decrease the area of influence 
of the grain boundaries. Besides, the segregation of cal- 
cium is strongly determined by the strain field. More 
calcium segregation is expected in tensile stress regions. 
Since the calcium is doped at yttrium sites, it has one 
negative charge. Therefore, the segregation of calcium 



causes extra negative defects to approach grain bound- 
aries such that the potential decreases quite significantly. 
Consequently, the hole content near grain boundaries in- 
creases, thus helping to screening out the effect of disor- 
der at the boundary. In summary, the nature of space 
charge is closely related to the oxygen content. Doping 
calcium at small angle tilt grain boundaries of YBCO 
leads to the increase of J c through increasing the hole 
content, enhancing negative potential regions due to the 
segregation of calcium, and helping to passivate disorder 
at the boundary. 

Note Added — After submission of this manuscript, 
ref . |39| reported the experimental measurement of po- 
tential curves near [001] tilt boundaries for both 
YBa 2 Cu307_(5 and Ca-doped YBa 2 Cu307_5. In partic- 
ular, the curves in Fig. 2 of ref. [33 were nicely predicted 
by Figs. 2 and 5 in this manuscript. (Indeed, the main 
body of this manuscript has been already completed as 
chapter 6 of ref.0. The Figs. 6.2 and 6.3 of ref.^ are 
in fact plotted here as Figs. 2 and 5.) 
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